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A search for the fJ-^fJ-^ and fJ-^fJ-^ decays is presented using ^ 300pb-i of pp 

' ^ ' colUsions at y(s = 7 TeV collected by the LHCb experiment at the Large Hadron Collider at 

D CERN. The measured upper limit for the branching ratio of the l^-^IJ- decay is ^{B^ — 

li^li^)< 1-3 (L6) X 10"** at 90 % (95 %) confidence level (CL), while in the case of the B^ 
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11+11- decay the measured upper Umit is .^(B^ ^+^")< 4.2 (5.1) x 10"^ at 90 % (95 %) CL. 
A combination with the 2010 dataset results in ^(B,. ^ ^+^")< 1.2 (1.5) x 10"^ at 90% (95 %) 
CL. 
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1. Introduction 

Measurements at low energies may provide interesting indirect constraints on the masses of 
particles that are too heavy to be produced directly. This is particularly true for Flavour Changing 
Neutral Currents (FCNC) processes which are highly suppressed in the Standard Model (SM) and 
can only occur through higher order diagrams. The SM prediction for the branching ratios (^) of 
the FCNC decays Bs IJ-^IJ-^ and — > M^M" have been computed to be (3.2 ±0.2) x 10"^ and 
(O.lOibO.Ol) X 10^^ respectively [1]. However New Physics (NP) contributions can significantly 
enhance these values. 

The best published limits from the Tevatron at 95% CL are obtained using 6.1 fb"^ by the DO 
collaboration [2], and using 2 fb^^ by the CDF collaboration [3]. The CDF collaboration has also 
presented a preliminary result [4] with 6.9 fb^' in which an excess of — )■ /i+/x^candidates is 
reported, compatible with a 1^^}^^)= (l-Slog) ^ ^0"^. 

The LHCb collaboration has previously obtained the limits SS{Bs lx^pL^)< 5.4 x 10^^ and 
^{Bd IX^IX^)< 1.5 X 10"^ at 95% CL based on 37 pb"' of luminosity collected in the 2010 run 
[5]. We present here a measurement based on 300 pb^' of integrated luminosity collected between 
March and June 20 11. 

2. Analysis strategy 

The general structure of the analysis is similar to the one described in Ref. [5] and is detailed 
in Ref. [6]. 

The selection procedure treats signal and control/normalization channels in the same way in 
order to minimize the systematic uncertainties. Assuming the SM branching ratio and the bb cross- 
section, measured within the LHCb acceptance, of a^^ = 75 it 14/xb [7], approximately 3.4 (0.32) 
B^ ^ H^IJ.^ (B^ — > jLt+/i^) events are expected to be reconstructed and selected in the analysed 
sample. 

After the selection, each event is given a probability to be signal or background in a two- 
dimensional space defined by two independent likelihoods: the invariant mass and the output of a 
Boosted Decision Tree (BDT) from the TMVA package [8]. The combination of variables entering 
the BDT is optimized using Monte Carlo (MC) simulation. The following variables have been 
used: the B lifetime, impact parameter and transverse momentum of the B, the minimum impact 
parameter significance of the muons, the distance of closest approach between the two muons, the 
degree of isolation of the two muons with respect to any other track in the event, the cosine of the 
polarization angle, the B isolation and the minimum pj of the muons. The BDT distribution is then 
transformed in order to be flat for the signal and peaked at for the background. 

The calibration of the invariant mass and the BDT likelihoods are obtained from data using 
control samples. The signal BDT shape is obtained from 5°^^ — )• events free from trigger 

biases while the background shape is obtained using sideband B^^^^ — jLt^/i^ candidates. The re- 
sulting distributions are shown in Fig. 1 . 

The parameters describing the invariant mass line shape of the signal are extracted from data 
using control samples. The average mass values are obtained from S° — )• K^n^ and — )• K^K 
exclusive samples. The and 5" mass resolutions are estimated by interpolating the ones obtained 
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Figure 1: BDT calibration for signal and background. 

with the dimuon resonances ^{2S) and T( 15, 25,35)) and cross-checked via a fit to the in- 
variant mass distribution of the Bj'^j — h^h^ inclusive decays and of the — > K^n^ exclusive 
decay. The interpolation yields g{B) = 24.6 ± 0.2stat ± 1 -Osyst- 

The number of expected signal events is obtained by normalizing to channels of known branch- 
ing ratios, B+ — J/yK^ , — >■ J/wP > ~^ K^n^ , that are selected in a way as similar as 
possible to the signal. 

The probability for a background event to have a given BDT and invariant mass value is ob- 
tained by a fit of the mass distribution of events in the mass sidebands, in bins of BDT. Different fit 
functions and mass ranges are used to compute the systematics uncertainties. The two-dimensional 
space formed by the invariant mass and BDT is binned, and for each bin we compute how many 
events are observed in data, how many signal events are expected for a given 3S hypothesis and lu- 
minosity, and how many background events are expected for a given luminosity. The compatibility 
of the observed distribution of events in all bins with the expected one for a given 3S hypothe- 
sis is computed using the CLs method [9], which allows to exclude a given hypothesis at a given 
confidence level. 

In order to avoid unconscious biases, the data in the mass region defined by M^o — 60MeV/c^ 
and Mgo + 60MeV/c^ have been blinded until the completion of the analysis. 

3. Results 

The distribution of events in the invariant mass versus BDT plane is reported Fig. 2. The 
expected limit at 90 (95) % CL for the S, iU+/i" is 0.8 (1.0) x 10"^ in the case of background 
only hypothesis. When adding signal events according to the SM branching fraction, these limits 
become 1.2 (1.5) x 10"^. The observed values for the jU+/i"channel is 1.3 (1.6) x 10"^ with 
a CLb value of 0.80. The observed events are in good agreement with the background expectations 
and the presence of — )• /^^/i^events according to SM predictions. 

For the Bd l^^l^^, the expected limit at 90 (95)% CL is 2.4 (3.1) x 10"^ in the case of 
background only hypothesis. The observed values is 4.2 (5.2) x 10^^ with a CLb value of 0.79. 
The comparison of the observed distribution of events with the expected background distribution 
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Figure 2: Distribution of selected dimuon events in the invariant mass vs BDT plane. The orange short- 
dashed (green long-dashed) lines indicate the ±60MeV/c^ search window around the (B"). 

results in a p-value of 20 % (21 %) for the — > M^M" (Bd — > M^M") decays. In the case 

of Bd — > M^M^, the slightly low p-value is due to an excess of the observed events in the most 
sensitive BDT bin with respect to the background expectations. A larger data sample will allow to 
clarify the situation. In the case of B^ — pL^pL^ , when a signal is included at the level expected in 
the Standard Model, the p-value increases to 50 %. 

Finally, the B^ — )■ pL^pL^ limit is combined with the one published from the 2010 data to obtain 
^(B,s M^M") < 1-2 (1.5) X 10"^ at 90 % (95 %) CL. This 90 % CL upper limit is still 3.8 times 
above the standard model prediction. 
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